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O ; Abstract 

Homogeneous nucleation of the new phase of one transition near a second phase transition is considered. The system 



(N 



^ | has two phase transitions, we study the nucleation of the new phase of one of these transitions under conditions such 



o 



that we are near or at the second phase transition. The second transition is an Ising-like transition and lies within 
the coexistence region of the first transition. It effects the formation of the new phase in two ways. The first is by 

■ reducing the nucleation barrier to direct nucleation. The second is by the system undergoing the second transition and 
CNJ \ transforming to a state in which the barrier to nucleation is greatly reduced. The second way occurs when the barrier to 

■ undergoing the second phase transition is less than that of the first phase transition, and is in accordance with Ostwald's 
rule. 
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1 Introduction 



| The formation of a new phase at a first-order phase transition is an activated process. A nucleus of the new phase must 
■ overcome a free energy barrier in order to form and then grow into the new phase. The rate at which such microscopic 
nuclei form is proportional to exp(—AF*/kT), where AF* is the height of the free energy barrier which the nuclei must 
d . overcome El . Here we consider nucleation of the new phase at a first-order phase transition and calculate AF* when 
Ph the system is near a second phase transition. One of the examples of nucleation near a second transition is nucleation 
I . of a crystalline phase at a first-order fluid-crystal transition near a fluid-fluid phase transition. Such a process has 
been observed in globular proteins whose phase diagrams show a metastable fluid-fluid transition within a strongly 
q first-order fluid-crystal transition (^,|^|. Crystallisation of globular proteins is the subject of interest as protein crystals 
' are needed in order to study their structure by X-ray diffraction ||-|6|] . It has been seen that globular proteins crystallize 
at temperatures near where we expect a metastable fluid-fluid critical point Numerical work on nucleation near 

a critical point has been done by Talanquer and Oxtoby 0, where they obtained nuclei with very large numbers of 
^ . molecules near the critical point. This followed pioneering computer simulations of ten Wolde and Frenkel G3], who 
' found an anomalously low AF* for the nucleation of a crystalline phase near the critical point of a metastable fluid- 
fluid transition. Other recent theoretical work on nucleation near a metastable transition may be found in Jll] - [l6| and 
' references therein. Earlier theoretical work by one of us |p^|l6[| showed that as the critical point is approached, the 
I , derivatives of the free energy barrier to nucleation with respect to the chemical potential and temperature diverge. 
y—i 1 Therefore the presence of the critical point causes a rapid drop in the free energy barrier to nucleation and so facilitates 
. nucleation. 

The metastable transition we consider here is an Ising-like or a vapour-liquid-like transition. An Ising-like transition 
is a transition from a phase with a negative magnetisation to a phase with a positive magnetisation. A vapour-liquid- 
£h like transition is a transition between two fluid phases differing in density. The metastable transition lies within the 
. coexistence region of the equilibrium transition. Thus when the system is in the coexistence region of the metastable 
<— | ' phase transition, it also lies within the coexistence region of the equilibrium transition. Then one of two new phases can 
O ■ nucleate: the stable phase can nucleate, as can the liquid phase or the analogue of the liquid phase, from the vapour 
' phase or the analogue of the vapour phase. Which of these occurs first depends on which free energy barrier is lower. 
^ . Thus the stable phase can be formed via two processes. In one process the free energy barrier to nucleation of the stable 
phase becomes small enough to allow nucleation, but the barrier to formation of the liquid phase is still large, then the 
stable phase nucleates. When the sizes of the barriers are reversed, then the formation of the stable phase occurs via 
C$ ' two steps: first the liquid phase nucleates, and second the stable phase nucleates in this liquid phase. 

Our finding that nucleation of a metastable transition can occur before nucleation of the stable phase, is just that 
observed by Ostwald more than 100 years ago jl7],[l8|. Ostwald's rule is that the phase that nucleates first is not 
necessarily the most thermodynamically stable phase, but is the one with closest free energy to the fluid phase. Later 
Stranski and Totomanow (l^] improved this rule and suggested that the nucleated phase is the one with the lowest free 
energy barrier to nucleation, and not the phase that is globally stable. 

Nucleation can be heterogeneous or homogeneous. In most circumstances suspended and dissolved impurities, as 
well as the solid boundaries provide sites for the formation of the new phase. This process is known as heterogeneous 
nucleation. However, in the absence of impurities or solid surfaces, small nuclei of the new phase are formed within the 
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bulk of the system. This is a homogeneous nucleation and this is what we are studying here. A free energy barrier must 
be overcome in order to form nuclei of a critical size, beyond which the new phase grows spontaneously fj],[20|. The 
rate at which critical nuclei of the new phase are formed is very sensitive to the height of the free energy barrier. The 



so-called classical nucleation theory was originated with the work of Volmer and Weber 20,2lJ| and originally developed 
for droplet condensation from supercooled vapours. The theory calculates the rate at which nuclei grow to a critical 
size, proportional to exp(— AF//fcT), where AF/ is the minimum work needed to form the critical nucleus. 

In this paper we carry out numerical calculations to study the homogeneous nucleation of the new phase at a first- 
order transition near or at the second phase transition. The second transition is an Ising-like transition from the phase 
with a negative magnetisation to the phase with a positive magnetisation. In this work we use magnetic language, in 
fluid language the transition is from a low density phase to a high density phase. We present the results of numerical 
calculations for the temperatures below the critical point. The calculations are for the excess order parameter Am* and 
the free energy barrier to nucleation of the new stable phase, AF*, in the critical nucleus. In the next section we study 
the theory of two types of homogeneous nucleation: nucleation of the new stable phase, and nucleation of the positive 
magnetisation phase from the negative magnetisation phase. We then discuss some of the numerical techniques we used 
in the calculations in section ||. In section || the results are presented and the last section is a conclusion. 

2 Theory 

2.1 Nucleation of a stable phase 

We study nucleation of a new phase at a first-order transition near a second transition. The second transition is an 
Ising-like transition with an order parameter m. We work below the critical or Curie temperature and so we have a 
transition from negative m to positive m on increasing the field h. h is as usual the field which couples to m. Following 
earlier work by one of us |l4| we split the nucleus of the new phase into two parts: the core and the fringe. That part 
of nucleus less than r c from the origin is the core and the part farther than r c is the fringe of the nucleus. We assume 
that the fringe of the nucleus is spherically symmetric. Therefore the order parameter m(r) is a function only of r, 
the distance from the centre of the nucleus. As the universal behaviour of the nucleus is derived from the fringe, we 
therefore substitute the core by a boundary condition on m(r) of the fringe. Thus for simplicity we fix m{r < r c ) = m c , 
where the subscript c stands for the core of the nucleus. m c is a positive constant independent of temperature and of 
h, and is essentially the value of the order parameter m in the nucleating phase. We require the excess free energy of 
the nucleus, AF, which is the free energy with a nucleus present minus that without a nucleus. Here we use a standard 
square-gradient functional for the fringe of the nucleus as a functional of the order parameter profile of the nucleus, 



AF = AF r 



Aoj{m) + k (Vm) 2 dr, (1) 



where the first term is the contribution from the core of the nucleus, AF C = |7rrj? Acj core , and the integral term is the 
contribution from the fringe of the nucleus. In Eq. ([[]) and all other equations, our free energies are in units of kT; k 
and T, Boltzmann's constant and the temperature, respectively. Au> core is a constant and is the excess free energy per 
unit volume of the core. As the core is close to the bulk equilibrium phase, Auj core < because the free energy is lower 
in the equilibrium phase. The gradient term in the integrand of Eq. (|l|) is due to the variations in space of m(r). This 
is the lowest order term in a gradient expansion and is adequate when m is slowly varying. The coefficient n is assumed 
to be a constant. In Eq. (|l]) Auj in the integrand is given by 

Auj = f(m) - f(m b ) - h(m - m b ), (2) 

where m;, is the order parameter in the bulk. Alu is the work required per unit volume to change the order parameter 
from its bulk value nib to m in the presence of an external field h. f(m) is the bulk Hclmholtz free energy per unit 
volume as a function of the order parameter m(r). The critical nucleus is at the top of the free energy barrier and 
therefore at the maximum of AF. So we set the functional derivative of AF with respect to the order parameter profile 
m(r), to zero, 

3Aw(m) n 2 



dm{r) 



2nV z m{r) = 0, r > r c , (3) 
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The total excess of the order parameter due to the presence of the nucleus, which is the total order parameter with the 
nucleus present minus that without the nucleus, is obtained as 

Am — J (m(r) — m^dr. (4) 
We employ Landau theory [^5],|26| for the bulk free energy, 

f(m,T) — —atm 2 + 6m 4 , (5) 

where a and b are constants, and t = T — T c with T c the critical temperature p6| . The bulk order parameter TOf,(/i,T) 
is obtained by solving 



dm 



(6) 



where co(m) — f(m) — hm. Then using Eqs. (^|J^) gives 



atm + 4bm 3 -h = 0. (7) 

When h = the bulk order parameter is that at coexistence and we call it m co . For h — 0, Eq. (Q) has just one real 
solution m co = for T > T c . For h — and T < T c , Eq. (0) has three solutions, and consequently the free energy 
has two minima with the same free energy at m co = ±| (—at/b) 1 ^ 2 , and a maximum at m — 0. In the presence of 
an external field h, Eq. (j?]) has one real solution for T > T c , which corresponds to a single minimum for uj{m). When 
T < T c and h ^ with \h\ not too large, there are two inequivalent minima in u>(m). For a fixed h, the absolute 
minimum corresponds to the equilibrium state and the higher minimum is a metastable state. 
To get the order parameter profile we substitute Eq. (^) in Eq. (||), which gives 



dr 2 r dr 2k 
for r > r c , subject to boundary conditions 



d 2 m 2 dm 1 , „ , , , . . 

(46m 3 + atm - h) = 0, (8) 



m(r 



(9) 



and 



m(r —> oo) = mj. (10) 

The former is the order parameter at the boundary between the fringe and the core, and the latter is the obvious 
condition that m tends towards its bulk value far from the nucleus. Eq. (||) is a nonlinear second-order differential 
equation which can not be solved analytically. We therefore use numerical methods to solve Eq. (j^) for the order 
parameter profile, and then substitute it in Eq. ([!]) to obtain the excess free energy of the nucleus. The barrier to 
nucleation is at the maximum of the excess free energy and so will occur at a radius r* given by dAF/dr c = 0. Thus we 
need to use a maximization method numerically to find r* and hence AF* , which is the free energy barrier to nucleation 
for a critical nucleus. This also gives m*(r) which is the order parameter profile of the critical nucleus. Substituting 
this profile in Eq. (^) we then find the excess order parameter of the critical nucleus Am* . 

2.2 Nucleation of +m phase from — m phase, Classical Nucleation Theory 

Starting in the — m phase and increasing m (or equivalently h), we cross the phase transition from the — m to the +m 
phase. Then our system can not only nucleate the equilibrium phase but also the +m phase. Then we are within both 
the — to to +to coexistence region and the coexistence region between the equilibrium phases. The nucleus of the +to 
phase must overcome a free energy barrier in order to form, and then it grows to the new phase. The free energy barrier 
to form such critical nucleus is given by 
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(11) 



which is the standard classical nucleation theory expression |1|,|2C|,|27|. This expression is valid for nucleation not too 
close to the spinodal [^8|. The subscript / indicates that the transition is an Ising-like phase transition. In Eq. ( |TT| ) a 
is the surface free energy of the equilibrium interface between +m and —to phases, and is given by 



a = 2 ^nAf(m)dm, (12) 



where A/(m) = /(to) — f(rrib) is the excess bulk free energy, Eq. (|[). 

It is worth mentioning that in this paper we have used a magnetic rather than a fluid language to describe the phase 
transitions. To switch from the magnetic language with a transition from — to to +to phase, to the fluid language with 
a vapour- liquid- like transition, the external field h becomes the chemical potential p minus that at the coexistence p co , 
the order parameter m becomes density p, the bulk order parameters —mb and +nib become the vapour and liquid 
densities, p v and pi, respectively, and the axis to = becomes the critical density p cp = (p v + pi)/2 in the coexistence 
curve. Also in the fluid language the total order parameter Am* is equivalent to the total number of molecules in the 
critical nucleus. 

3 Numerical Techniques 

In the fringe far from the core, m(r) is near the bulk order parameter mb, we therefore use a Taylor expansion of Au(m) 
about to = mb, 

Aw(m) = -X _1 (m-m&) 2 + --- (13) 

-^-l= x -\m-m b ) + --- , (14) 
where Au> and its first derivative are zero at m = mb- Here \ is t ne response function of m, defined by 



X 



i _ d 2 f 



dm 2 

where /(to) is the bulk free energy Eq. (||). Substituting Eq. (14) in Eq. (||) gives 



(15) 



X 1 {m — mb) — 2kV 2 (to — mj) — 0, (16) 

which has a solution of Yukawa form 

m(r)=m b + B C -^tlM i (17) 

r 

for large values of r. To obtain Eq. jl7|) we used the boundary condition Eq. ([Tol). In Eq. (|l7| ) B is a constant value 
which we will fix later, and £ is the correlation length for m(r) which is given by 

£=(2 KX ) 1/2 . (18) 

The solution of m(r) for large values of r has the form of Eq. (|l7|). In numerical calculations we therefore use a large 
distance from the center of the core, R, to be the outer boundary of our numerical integration of Eq. d|), and use 
Eq. @ beyond R. 

For our numerical calculations we employ the Runge-Kutta method to solve the second-order differential Eq. (||) 
for the order parameter profile, subject to boundary conditions at r c and R, Eq. (||) and Eq. (p~7|) , respectively. The 
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constant value B in Eq. ( ^7[) is then obtained by finding the root of the nonlinear equation m(r c ) — m c = 0, for which 
we used Ridders' method [£9|. Substituting the order parameter profile m(r) in Eq. (Q) gives the excess free energy of 
the nucleus. The barrier to nucleation is at a maximum of AF and so will occur at a radius r* given by dAF/dr c = 0. 
We then use a maximization method to obtain r* and consequently to* (r) and the barrier to nucleation AF* . 

In the numerical calculations we used a large radius R for the boundary condition of the fringe of the nucleus. We 
need such a value to enable us to carry out the numerical calculations. We calculate the contribution of distances larger 
than R analytically and add it to the corresponding numerical quantity for the range r c to R, as a correction term. The 
correction term for the excess order parameter of the critical nucleus, Am*, can be calculated using Eq. (Q) with the 
limits for the integral from R to oo. Then by substituting Eq. ([l7]) for m(r) for all values of r > R, the correction term 
is 

Am*, = 4itkB£,(R + £)e- R K, (19) 

where ct stands for the correction term. The correction term for the free energy barrier to nucleation AF* , is calculated 
using the integral term in Eq. (Q) for all values of r from R to oo. Then by substituting Eq. ( |T^ ) in Au>(m) and using 
the normal gradient of Eq. (^) for Vm(r), the correction term is given by 

AF* = 4™B 2 (i + e" 2 ^, (20) 

which we have used the definition of £, Eq. (p"8|). 

The solution of Eq. ( p"7| ) for the order parameter for r > R is derived using the Taylor expansion of Aw(m) about 
to = nib, Eq. (|l^). In this expansion we have used up to the second-order term of (m — mb) and have ignored the higher 
orders. We therefore need to obtain the conditions under which we are allowed to use this approximation. To do so we 
see that in Eq. ( |l3| ) by ignoring terms of orders greater than 2, we have assumed that these values are less than the 
second-order term. The biggest of all ignored terms is the third order one, we have therefore assumed that 



9 3 Aw(m) 



m — nib) 1 _i 



3! 2 

m—m^ 



< -X' L (rn-m b y. (21) 



dm 3 

Using Eqs. (§,|lj|) in Eq. @ gives 

— m b ?(m(R)-m b )<*:l, (22) 

K 

as the requirement that must be satisfied. Here m(R) is given by Eq. ( |l7| ) at r = R. 

4 Results 

The free energy barriers to nucleation of the stable phase and of the +m phase are AF* and AF/ , respectively. In each 
nucleation the formation of the new phase occurs when the free energy barrier is of roughly 30fcT. The figure 30kT 
is rather arbitrary [0 but the nucleation rate is of order, number of molecules times exp(— AF*/kT) divided by the 
characteristic time scale of the solution. The characteristic time scale of protein solutions is of order 1/is ||. Thus a 
volume of the solution containing 10 9 protein molecules has a nucleation rate of order 10 9 x e -30 x 10 6 ~ lOs^ 1 . Thus 
in a volume of the solution containing a billion molecules, nuclei will appear in a fraction of a second. When the barrier 
reaches 30fcT the — m phase is unstable and the stable +m phase nucleates immediately. A barrier of about twice this, 
60/cT, is required for the —to phase to be stable over long times, i.e., for it to be metastable. Homogeneous nucleation 
of a new stable phase near a second phase transition has been studied before |p|-|l6[|. Here we carried out numerical 
calculations for Am* and AF*. Near the critical point the scaling of the size of the critical nucleus is [14-[lq], 



Am* ~ \tr 9± (h/\tf s ) , (23) 

and that of the free energy is 

AF S * ~ \tff ± {h/\tf s ) , (24) 
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Figure 1: Am* versus h at two temperatures below the critical temperature, when k = 1 and Auj core = —0.5. Solid and 
long-dashed curves correspond to t = —0.5 and t — —0.25, respectively. 



for both temperatures above and below the critical temperature. Here ± are related to t > and t < 0, 7, j3 and 8 are 
the usual critical exponents, and take values of 1, 1/2 and 3 within mean-field theory. In Eq. ( p4| ) subscript s stands for 
singular part of the fringe. These results were predicted in earlier work |l4|-|l6| and we confirmed them by numerical 
calculations. As the critical point is approached, the magnetisation of the critical nucleus Am* diverges, and so does 
the derivative of the free energy barrier to nucleation AF*, with respect to the temperature and to h. We are also 
interested in the behaviour of Am* and AF* far from the critical point. 

4.1 Far from the critical point 

Here we study the total magnetisation and the free energy of the critical nucleus, Am* and AF*, respectively, far from 
the critical point — \t\ not small. In all the calculations we have fixed a = b = m c = 1. In Figs. |l| and || there are graphs 
of Am* and AF*, respectively, for two different temperatures below and far from the critical temperature, t = —0.25 
and t = —0.5 with k = 1 and Auj core — —0.5. For each temperature, when h < the system is in a state with a negative 
order parameter m < less than that at the coexistence. When h > the system can be found in a state with a positive 
order parameter m > greater than that at the coexistence. For h > there are also states with m < between the 
— m branch of the coexistence curve and the spinodal. These states are within two coexistence regions: the coexistence 
region between — m and +m phases, and that between the equilibrium phases. Therefore in both Figs. [I] and H for each 
temperature there are values of h > for which Am* and AF* are double valued. These values correspond to the — m 
phase between the coexistence and spinodal curves, and the +m phase with a lower Am* and AF*. 

In Fig. [I] for each temperature as h increases the curves on the right-hand side decrease. Those on the left-hand 
side decrease first and then increase rapidly. Decreasing Am* is due to the fact that for a higher h, the bulk order 
parameter m\, is higher and thus it leads to a lower excess order parameter Am* in Eq. (||). However the large and 
rapidly increasing values of Am* are due to the proximity to the spinodal. As we approach the spinodal the response 
function of the order parameter, \ defined as in Eq. (|I1f ) diverges. Therefore the total excess order parameter Am* 
which scales as \ |l4|-|l(| will diverge as well. As the temperature approaches T c (\t\ decreases), Am* diverges at a lower 
value of h. This is because, as T approaches T c , the gap between the coexistence and the spinodal curves decreases. In 
Fig. U for each temperature as h increases AF* decreases, due to the increase in the bulk order parameter m&. The phase 
with the higher order parameter must overcome a lower free energy barrier AF* to form a nucleus of the equilibrium 
phase, whose order parameter we have set equal to 1. Also in Fig. @, AF* varies rapidly as it approaches the spinodal. 
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Figure 2: AF* versus h for two different temperatures, when k = 1 and Acj core = —0.5. Solid and long-dashed curves 
correspond to t = —0.5 and t = —0.25, respectively. 



This is due to what is called the nucleation theorem pfj|-|32[] and states that the derivative of the free energy barrier 
with respect to h is equal to minus the excess m, dAF*/dh = —Am* It shows that the larger Am* is, the 

more rapidly the free energy barrier varies with h. As the spinodal is approached Am* diverges and thus AF* drops 
rapidly with h. When a transition from — m to +m phase occurs both Am* and AF* in Figs. [I] and ^| jump to lower 
values, because of the increase in m^. The height of this reduction is temperature-dependent because the bulk order 
parameters m co = ±y/ (— 1/4) are functions of t = T — T c . As the temperature decreases, the difference between these 
two order parameters increases, and therefore the height of the drop in both Am* and AF* increases. 

For each temperature as h increases AF/ decreases, Eq. (pT[), and for energy barriers of about 30fcT a transition of 
— m phase to +m phase with lower values of Am* and AF* occurs. The cross marks on the graphs in Figs, [j] and |^ are 
the points where AF/ = 30fcT. Also in the region with h > and m < the solutions to Am* and AF* disappear at 
a value of h before the spinodal. The circle marks on the graphs are the points for which the solutions disappear. For 
this value of h in our numerical calculations the solution to the order parameter profile of the critical nucleus, m*(r), 
disappears. The solution to m*(r) is obtained by solving Eq. ([?[) subject to the boundary condition m*(r = r*) = m c . 
In the region between the — m branch of the coexistence and the spinodal with h > 0, Eq. (^) gives two solutions for 
m*(r), out of which we choose the one with lower free energy barrier, AF*. As h increases, these two solutions approach 
each other until they become equal at a specific value of h > 0. For any value of h greater than this, then there is no 
solution. This means that the critical nucleus does not exist anymore and there is no well defined free energy barrier 
to nucleation. As we see in Figs, [j] and ^| for each temperature the points with AF/ = 30fcF (cross marks) are at the 
left side of those for which the solutions disappear (circle marks) . As for energy barriers of AF/ = 30kT a transition of 
— m to +m phase occurs, therefore for these parameter values as the system approaches the cross marks, the —in phase 
will transform to +m phase. Thus in the experiments the states after the cross marks are not accessible. As we see in 
Figs. and || as the temperature approaches T c (\t\ decreases), the solution disappears at a lower value of h > 0. This 
is because as T approaches T c , the coexistence and the spinodal curves approach each other. 

The stable phase can be attained via two processes |3^]. In one process AF* drops below 30fcT when AF/ is still 
large, then the stable phase is reached in one step through nucleating in — m phase. In the second process AF/ drops 
below 30fcT when AF* is still large, then the stable phase will be reached in two steps: first the +m phase nucleates 
in the — m phase and grows and then second, the stable phase nucleates in this +m phase. The fact that the nucleated 
phase is the one which has lower free energy barrier and not necessarily the stable phase, is in accordance with Ostwald's 
step rule. In Fig. ^ with Auj core = —0.5 the barrier AF* is always greater than 30fcT when AF/ drops below 30fcT. 
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Figure 3: AF* versus h for two Auj core , when k — 1 and t — —0.25. Solid and long-dashed curves correspond to 
Auj core = —0.75 and Auj core = —1.25, respectively. 



Thus on increasing h the +m phase nucleates before the equilibrium phase. The equilibrium phase then nucleates in the 
+m phase. But as Auj core < 0, when Auj core becomes more negative then AF* becomes more negative and therefore 
AF* in Eq. (jl|) drops. So for sufficiently negative values of Auj core , AF* drops below 30fcT before AF/ does and 
the stable phase nucleates before the +m phase. To show this, in Fig. || graphs of AF* are plotted for two values of 
Auj core — —0.75 and Aoj core = —1.25, when t = —0.25. Just as in the previous figures, the cross marks on each graph 
correspond to where AF/ = 30kT and the circle marks are where the solutions disappear. For Aoj core — —0.75, AF* 
is greater than 30fcT when AF/ drops below 30fcT. Thus the +m phase nucleates first and then the equilibrium phase 
nucleates in the +m phase. For Au) core = —1.25 however AF* drops below 30fcT before AF/ does. Thus the equilibrium 
phase nucleates before the +m phase. The values of h for which AF/ = 30fcT (the cross marks) are independent of 
Acu core and are the same on both graphs. This is because AF/ in Eqs. ( |ll|Jl2] ) depends only on the temperature and 
on h, and is independent of Aui core . 

The two step route to the equilibrium phase that we find is not dissimilar to that proposed for the crystallisation 
of some polymer melts by Olmsted et. al. p4| . They study a system with the same phase diagram as considered here 
but their equivalent of the — m to +m phase transition proceeds via spinodal decomposition and not nucleation as here. 
Thus the results here apply to liquids and colloidal suspensions, where the new phase nucleates, whereas Olmsted et. al.'s 
approach is more applicable to polymers. 

We are also interested to study where the free energy barrier AF/ = 30fcT, and a nucleation from the — m phase to 
the +77i phase occurs. Up to now we considered systems with k = 1. We also like to know the effect of increasing k on 
the calculations. In Figs. || and there are diagrams of t versus m, and t versus h, respectively for three different values 
of k when Aui core — —0.5. In Fig. |i| the solid (outermost) and long-dashed (innermost) curves represent the coexistence 
and the spinodal, respectively. The other curves in Fig. |] between the coexistence and spinodal curves correspond to 
the points with energy barriers of AF/ = 30fcT for different k. From left to right, dotted, dashed and dot-dashed curves 
are for n — 1, 1.5 and 2, respectively. In Fig. ^| the solid line h = along axis t corresponds to coexistence, and the 
long-dashed curve (the outermost) represents spinodal. The other curves in Fig. || between the coexistence and the 
spinodal correspond to the points with AF/ = 30kT for different n. From left to right, dotted, dashed and dot-dashed 
curves represents k = 1, 1.5 and 2, respectively. As k increases there is a shift towards larger values of both m and 
h for AFj — 30kT. This shift can be explained as follows. For a fixed temperature and external field h, a higher k 
yields a higher free energy barrier to nucleation of the +m phase AF/ ([Tl| , [l2] ). Now for each fixed temperature and k, 
as h increases the energy barrier AF/, Eq. ([ll]), decreases until it reaches energies of about 30fcT. Therefore at each 
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Figure 4: Phase diagram for t versus m. Solid and long-dashed lines are coexistence and spinodal, respectively. The 
others correspond to AF 7 * — 30fcT, dotted, dashed and dot-dashed lines are for k — 1, 1.5 and 2, respectively. 



fixed temperature, for a higher k the nucleation occurs at a higher external field h and consequently at a higher order 
parameter m. Increasing k also increases AF* Eq. ([!]), thus increasing n pushes nucleation to both the +m and the 
stable phase to larger values of h. 

5 Conclusions 

Homogeneous nucleation of a new phase at a first-order transition near a second transition is considered. The second 
transition is an Ising-like or a vapour-liquid-like transition, and lies within the coexistence region of the equilibrium 
transition. Numerical calculations were performed to calculate the excess order parameter Am* and the free energy 
barrier to nucleation AF*, of the critical nucleus. For the states in the negative magnetisation phase with h > the 
system is within the coexistence regions of both transitions. Then in these states either the positive magnetisation 
phase or the equilibrium phase can nucleate. The free energy barriers to nucleation are equal to AF/ and AF*, for the 
former and the latter, respectively. The phase which nucleates is the one with the lower free energy barrier. This is in 
accordance with Ostwald's step rule, which states that the nucleated phase is not necessarily the stable phase, but is 
the one with the lowest free energy barrier to nucleation. The stable phase can be attained via two processes. In the 
first process AF* becomes small enough to allow nucleation when AF 7 * is still large, then the stable phase nucleates. In 
the second process AF/ becomes small enough to allow nucleation when AF* is still large, then the stable phase will 
be reached in two steps: first the positive magnetisation phase nucleates and grows and then second, the stable phase 
nucleates in this positive magnetisation phase. 

When the Ising-like phase transition occurs, both Am* and AF* jump to lower values, due to the increase in the 
bulk order parameter. As the bulk order parameter at coexistence is a function of the temperature, then the height of 
the jump is temperature-dependent. For a lower temperature, the jump in the order parameter increases and therefore 
the jumps in Am* and the barrier to nucleation AF* increase. This may mean a jump from a nucleation barrier too 
high to permit significant nucleation, to a very small barrier that results in rapid nucleation with many nuclei forming. 
This effect may cause problems for crystallographers wanting to achieve a small but non-zero rate of nucleation in order 
to obtain a few large crystals. We also observed that as the spinodal is approached, Am* diverges and the free energy 
barrier AF* varies rapidly with h. 
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Figure 5: Phase diagram for t versus h. Coexistence lies along the t axis. Long-dashed curve is the spinodal. The others 
correspond to AFf = 30fcT, dotted, dashed and dot-dashed lines are for k = 1, 1.5 and 2, respectively. 
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